The nature of self-localization in optical lattices with local dissipation 
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We analyze the nature of a novel type of self-trapping transition called self-localization (SL) of 
Bose-Einstein condensates in one-dimensional optical lattices in the presence of weak local dissipa- 
tion. SL has recently been observed in several studies based upon the discrete nonlinear Schrodinger 
equation (DNLS), however, its origin is hitherto an open question. We show that SL is based upon 
a self-trapping crossover in the system. Furthermore, we establish that the origin of the crossover 
is the Peierls-Nabarro barrier, an energy threshold describing the stability of self-trapped states. 
Beyond the mean-field description the crossover becomes even sharper which is also reflected by a 
sudden change of the coherence of the condensate. While the crossover can be readily studied in 
ongoing experiments in deep optical lattices, our results allow for the preparation of robust and 
long-time coherent quantum states. 
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Dissipation is typically known to represent a major ob- 
stacle in the coherent control of quantum systems. How- 
ever, in recent years, a strong interest in engineered dissi- 
pation has evolved, where dissipation has been used as a 
tool for quantum state preparation [TJ [5] as well as quan- 
tum information processing and entanglement generation 
[3j and to induce self-trapping (ST) [4-7_. Bose-Einstein 
condensates (BECs) have been shown to support a vari- 
ety of different kinds of ST, both in the continuous case 
(such as bright and dark solitons |MT5]) and in discrete 
systems |14fE^ . A particularly high level of control has 
been achieved in a two-mode BEC \^ , where ST can also 
be induced by local dissipation which can even repurify 
a BEC p]. 

A novel self-trapping transition coined 'self- 
localization' has been observed numerically in several 
studies based upon the DNLS in the presence of weak 
boundary dissipation in one-dimensional deep optical 
lattices |24ff?7] . In contrast to self-trapping, where a sys- 
tem is either prepared in a self-trapped state |T4t [TSl [28] 
or driven towards it ^ (6j ^ , SL is a mechanism where 
in presence of weak local or boundary dissipation a very 
general initially diffusive state leads to the formation 
of one or more discrete breathers (DBs, see |291 [30] 
for an overview). However, SL was only found, if the 
atomic interaction strength exceeds a critical value |26j . 
While the phenomenology of SL has been studied, the 
mechanisms that lead to this transition have remained 
unknown up to now. 

In this letter, we propose a mechanism for SL allowing 
us to give an explicit formula for an upper bound estimate 
of the SL threshold for the DNLS in excellent agreement 
with the numerical findings of [26J. The mechanism is 
based on a 'crossover' which surprisingly becomes much 
sharper when quantum corrections beyond the mean-field 
description are included, which is observed, e.g., in the 
condensate fraction of the system. 

To understand the nature of SL it is essential to note 
that the fixed point corresponding to the DB state into 
which the initial condition collapses does not undergo a 



bifurcation itself. On the contrary, using standard meth- 
ods |31l[51] the bright breather fixed point can easily be 
numerically found to exist and to be linearly stable for 
all positive nonlinearity strengths (and in the presence of 
boundary dissipation it becomes attractive). Linear sta- 
bility analysis therefore does not suffice to understand 
the SL transition. The underlying idea of our approach 
is that near the SL threshold a single strong, localized 
fluctuation of the number of atoms locally brings the sys- 
tem's state into the vicinity of the DB fixed point in phase 
space. In the simplest and most likely event this strong 
increase in the number of atoms happens on a single site 
that will become the center of the DB to be formed. 

We therefore study first, how a single site excitation 
can lead to the formation of a DB and find that there 
exists a distinct nonlinearity strength at which this ini- 
tial condition crosses over into a self-trapped state. We 
show that the origin of this ST crossover is an en- 
ergy threshold describing the stability of self-trapped 
states (called the Peierls-Nabarro (PN) energy barrier 
[221 [55]). Secondly, we statistically estimate the critical 
nonlinearity at the onset of SL by studying the proba- 
bility that a fluctuation in a diffusive state exceeds this 
ST crossover and leads to the formation of a breather. 
The ST crossover and SL should not only be observable 
for BECs but as well, e.g., in coupled nonlinear optical 
waveguides |29l [55] . 

Though SL is related to ST, it is distinguished by the 
way in which a stable (or metastable) and spatially lo- 
calized state is reached. There are several ways to ob- 
tain self-trapping of BECs in optical lattices which we 
classify into three types. Type I {'static preparation'): 
The quantum system is prepared in (or sufficiently close 
to) a self-trapped state. This has been realized in var- 
ious experiments \T\\ [M] ITSl [37] . Using a variational 
approach, a phase diagram has been calculated, that de- 
scribes the transition from diffusion to ST for an ini- 
tial Gaussian wave packet [1.9, 20., which, however, does 
not account for SL. Note that recent numerics for the 
DNLS P7I rather contradicts the phase diagram in [19]. 
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Type II ( 'dynamical preparation') : Another route to ST 
is to apply a strong local dissipation pulse, which can 
depopulate one or more sites and create a stable isolated 
peak or vacancy [H [51 [71 [53] (leading to the formation 
of a bright or dark breather). In particular, spatially re- 
solved dissipative manipulation in an optical lattice using 
an electron beam with single-site addressability has been 
demonstrated ^4^. Type III {'self -localization'): A third 
way to generate self-trapping is SL, where the system 
prepared in a random (generic) state in the presence of 
boundary or other local dissipation dynamically forms 
one or more DBs. Surprisingly, the locations where DBs 
form are not determined by the location of the leak [2$- 
[27] . In absence of boundary or local dissipation SL was 
not observed. 

Consider the Bose-Hubbard Hamiltonian in the mean- 
field description ^38, ,39j 

M J Af-l 

i/ = C/^|^„|^-- ^(CV'„+i+c.c.) (1) 

i—n n—1 

with on-site interaction U, tunneling rate J, lattice index 
n = 1 . . . M, where M denotes the number of lattice sites. 
Including boundary dissipation, the mean-field equations 
of motion are given by the dissipative DNLS (see [51 [71 [30] 
for a derivation of the loss term) 

ilpn = i|V-'n|^V'«-^(^/'n-l+V'n+l)-«7V'n('^«4+^n,M) (2) 

with Ti = 1, nonlinearity L = {2U/J)J\f, dissipation 
rate 7, total number of atoms J\f and the normalization 
IV'nP = 1- We introduce a measure of the local non- 
linearity V-""""^ = {2UM/J)Nn , where iV„ = \ipn\'^ is the 
relative number of atoms (also referred to as the norm) 
at site n. 

Let us first consider the dissipationless case with the 
following initial condition, where all atoms are located at 
site c, given by 

i^n{t = Q)=5nc- (3) 

In which range of the nonlinearity will the atomic popula- 
tion stay self-trapped at site c (resulting in the formation 
of a DB)? In Fig. [T|the evolution of the particle density 
is shown. For L = 1.6 (Fig. [ija)) the particle density 
initially decays exponentially in time and then populates 
the whole lattice evenly. In contrast, a completely differ- 
ent behavior is observed for L = 2.4 in Fig. [ijb), where 
the initial condition relaxes into a ST state which is ex- 
ponentially localized in space. A necessary condition for 
ST is LJ^'^^' > Leo, where Leo is the value of the nonlin- 
earity at the crossover that separates the diffusive from 
the ST regime. We define that ST is encountered, if 
min |'0c(t > r)p > a for large T, which is independent 
of T once a breather has formed. The value for a can be 
estimated via the position of a saddle point (in the so- 
called Peierls-Nabarro energy landscape, see below) that 
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Figure 1. A ST crossover for a (5-like initial condition, where 
all atoms are located at a single site c, is found both beyond 
and within the mean-field description (depicted here for the 
DNLS). The color code shows the normalized atomic density 
\^n{t)\^ . (a) Below the crossover, for L — 1.6 < Leo, the 
localized peak at t = decays exponentially fast, (b) Above 
the crossover (shown is L = 2.4 > Leo) a discrete breather 
forms. The particle density is stable and decays exponentially 
in space away from the center. About 85% of the atoms are 
located in three sites after time t = As. Other parameters are 
M = 101, c = 51, J = 10 Hz and 7 = 0. 



dictates the stability of the DB [22j, which is shown in 
Fig. [2|^a). In the limit L— )-oo the saddle point is found 
analytically at A2 = 1/2 (and Ni = N^^ = 1/4) [22J, we 
therefore estimate a — 1/2. Starting with initial condi- 
tion ([3| and choosing T = Is, the crossover from diffu- 
sion to ST is numerically found to be at LJ?^'" = 2.2463. 
Integration times were at least lOT. 

In the following, we will examine the ST crossover, for 
which we make use of a general concept called the FN en- 
ergy barrier. It is given by the energy difference \Ei,—Ee\, 
where L^f, is the total energy of a DB centered at a sin- 
gle lattice site and is the energy of a more extended 
breather centered between two lattice sites |3S] . The FN 
barrier is based on the notion, that due to continuity, 
the process of translating a localized object with energy 
Eb from one lattice site to the adjacent one involves an 
intermediate state with different energy E,,. 

We will connect the ST crossover to the stability of a 
DB. It has been shown that the stability of a DB can 
be well-described via a reduced problem of only few de- 
grees of freedom [29J, which reflects the fact that the 
breather is highly (exponentially) localized. This 'local 
Ansatz' has been further developed analytically in a local 
trimer (which is a subsystem consisting of three sites) on 
the so-called FN energy landscape [35], which is defined 
by i?pN = max50^^. (i/), with V'n = A/^exp(i(?i„) and 
54>ij = 4>i- 4>j [41 . The FN landscape reads [11] 

i/pN = ^ {Nl + + A|) + ( v/A^ + v^) . (4) 

Figure [2|^a) shows the FN landscape of the trimer at 
the ST crossover. The bright DB, which is linearly stable 
|42j, is located in the top 'eye' of the energy landscape. 
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The two saddle points just below N2 = 1/2 (related to 
a migration of the DB from site 2 to site 1 and 3 re- 
spectively) are connected to the PN barrier and the total 
energy threshold dictating the breather stability is given 
by 122j 

The energy of a bright breather Ei, is a maximum of the 
total energy E of the trimer. As long as the total energy 
of the local trimer Ep-^ < E < Ei, is above the threshold, 
a breather remains pinned to a lattice site. The total 
energy for the initial condition (|3| reads E{L) = L/2, 
which can be seen directly from Eq. (jlj as the energy is 
measured in units of the tunneling rate J (cf. Eq. (|2|). 
Hence, the crossover L^o is reached, when -BpN (Eq- ^) 
is equal to L/2, and we obtain 

Ll ~ 2Lt - Ll + Ll - + ? = . (6) 

We find Leo = 2.2469, in excellent agreement with the nu- 
merical value. This result means that the ST crossover, 
which is observed in a one-dimensional optical lattice, 
can be described with high degree of accuracy by the PN 
barrier of a local trimer. Given that the PN barrier de- 
scribes the stability of self-trapped states in a very broad 
context, we expect that the three different types to ob- 
tain ST (static, dynamical and self-localized) in discrete 
systems eventually are related to the PN barrier. 

The general behavior near the ST crossover is depicted 
in Fig. [2]jb). The PN barrier bends off the total energy 
line for increasing L > Leo, which leads to a growing area 
of stability (given by Eq{L) > Ep-^{L)). In the limit 
L — >■ (X), the initial total energy Eq (red line) asymp- 
totically approaches the total energy Ej, of the bright 
breather (blue thick line) The exact breather en- 

ergy Eb{L), here for M = 101 sites can be calculated 
numerically using standard methods (such as the anti- 
continuous limit [51^: 'ST, 13 ) , while we applied a different 
iterative approach [34 . 

Beyond mean- field.- To study, how the ST crossover 
manifests beyond the mean-field description, we use the 
Bogoliubov Backreaction (BBR) method, which includes 
higher-order correlation functions and enables a consis- 
tent calculation of the condensate fraction of the BEG 
|45|, I^J . The condensate fraction c/ is the fraction of the 
number of condensed atoms and is given by the largest 
eigenvalue of the reduced single-particle density matrix 
|38| . The BBR method has previously been generalized 
to the dissipative case [6] [7] . 

In Fig. ^c) the minimum remaining number of atoms 
(normalized to 1) at the central site are shown for M — 9 
sites and initial condition ^ using BBR and compared 
to the mean-field result. Boundary dissipation was ap- 
plied in both cases, reducing refiections from the edges 
of the lattice. By including quantum corrections, the ST 
crossover becomes much sharper which is also refiected 
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Figure 2. (a) The PN energy landscape -ffpN exhibits a 
bright DB (located in the top 'eye' on the PN landscape 
at N2 = 0.902 and N-^ ^ N3 ^ 0.049) and two degener- 
ate saddle points that mark the boundary of stability of the 
DB. The color code shows the PN landscape (Eq. Q) for 
L = Leo = 2.2469. (b) The ST crossover is found at the 
crossing (dashed line) of the total energy Eq = L/2 (red line) 
with the PN energy barrier -Epn (Eq. (|5|, green line). The 
corresponding energy at the crossing is the total energy of the 
saddle points shown in (a). For L — >■ 00, Eq asymptotically 
approaches the total energy of the bright breather Eb (blue 
line), (c-d) Including quantum corrections using the BBR 
method (blue line), the ST crossover exhibits a much sharper 
transition compared to the mean-field result (red line), here 
shown for M = 9 sites. We report the minimum number of 
atoms at the central site min(A'^5) and the condensate fraction 
min(c/(t)). The dashed line in (c-d) marks the crossover at 
Leo for an infinite lattice. Other parameters are 7 = 0.5 Hz, 
N = 200 atoms, J = 10 Hz. The minima were determined in 
the interval t £ [0.25, 0.5] sec. 



by a jump in the condensate fraction, see the blue curve 
in Fig. [2|d). In contrast, the mean-field dynamics based 
upon the DNLS per se assumes a pure BEG, i.e., c/ = 1 
(red line). While stable motion above the crossover al- 
lows for long-time coherence, unstable motion below the 
crossover leads to depletion of the condensate ^ 
profound understanding of the ST crossover therefore 
might be very useful for controlled quantum state prepa- 
ration using spatially localized initial conditions, such as 
Eq. (§. 

Self-localization. - We now turn our focus to SL, where 
the dynamics finds self-trapped states by itself in pres- 
ence of weak boundary dissipation [24H27j , resembling a 
phase transition |26J. To consistently investigate the dy- 
namics in different lattices sizes M, we require the initial 
density p = Af/M to be constant. Rescaling L accord- 
ingly results in an effective nonlinearity A — L/M [26J. 
Starting with a homogeneous initial condition with equal 
norm on all lattice sites and random phases, the tran- 
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Figure 3. The SL transition at A(, (red line) given by Eq. ( |11[ ) 
is an upper bound numerical results. The data (black circles) 
was extracted for boundary dissipation rate 7 = 0.2 Hz from 
Fig. 1 in |26| . where a sharp drop of the 'participation ratio' 
clearly marks the SL transition. The inset shows the proba- 
bility for detecting a norm x (Eq. ^) which exhibits a sharp 
transition and becomes a step function in the limit M 00. 



sition to SL has been observed at a critical interaction 
strength A;, for which we will derive an explicit expression 
in the following. The condition for self-trapping reads 
2^iocai ^ 2,iV„ = AMNn > Leo- The critical nonhnearity 
Af, for the dynamical formation of a breather is obtained 



for L 



local 



Leo J hence we find 



A, 



Leo 
MNr, 



(7) 



As the only unknown quantity in Eq. ([7| is the maxiniuni 
single site norm N^, calculating Af, is reduced to a very 
general question: What is the probability to find a site 
with norm larger than a given value N in the optical lat- 
tice? In the diffusive regime, the probability distribution 
of norms x in the lattice is w{x) = M exp{—Mx) [26J. 
The probability that the norm at a certain site is smaller 
than X is 



P{N <x)= w{x')dx' = 1 - e" 



(8) 



Assuming that the populations at the M sites are inde- 
pendent from each other, the probability that at least 
one site has a norm larger than x reads 



Pm[N > x) = 1 - [1 



-Mx}M 



(9) 



which approaches a step function for M — > 00 (see inset 
of Fig. |3|. Thus, the largest norm that is found in the 
diffusive regime is given for large M by Pm{N > x) w 
1/2 (red dashed line in Fig. |3|. Insertion into Eq. ([9| 
yields 



Nm = X = Inl 



]/M . 



(10) 



-1- (l/2)i/A/J 
With Eq.[7]the SL transition is found to be at the critical 



A, 



Lr 



(11) 



(1/2) 



which is shown in Fig. [s] (red line). As in deriving 
Eq. ( 10 1 it was assumed that the populations at the M 
sites are independent, we have effectively calculated an 
upper bound to A;,, in excellent agreement with the nu- 
merical results in [26] (shown as black circles in Fig.[3|. 

Experiments. - We expect that the ST crossover can be 
readily studied in current experiments. The initial con- 
dition ([3]) relates to a BEC cloud at a single lattice site, 
while the interatomic interaction can typically be con- 
trolled with high accuracy, e.g. via a Feshbach resonance. 
For large lattices (such that reflections from the edges are 
negligible), finite size effects won't play a role (see Fig.[lJ. 
In small lattices, refiections from the edges may slightly 
alter L^o but the system still exhibits the ST crossover 
both in the mean-field limit and beyond. In contrast, ob- 
serving SL in optical lattices is much more delicate. Weak 
local/boundary dissipation is a necessary prerequisite to 
observe SL, but could already be sufficiently present due 
to the inevitable coupling of the trapped condensate to 
the environment. The experiment, however, needs to al- 
low for sufficient propagation time so that ST can form, 
in the course of which chaotic dynamics and dynamical 
instabilities typically lead to depletion of the condensate 
|37[ l47l |48] . A remedy could be to reduce the relevant 
timescale by preparing an initial condition, that has more 
than exponential probability for high norms. SL repre- 
sents an alternative way to induce localization where the 
preparation of initial wave packets is not necessary. 

In conclusion, we analyzed the nature of SL in optical 
lattices, which previously has been observed phenomeno- 
logically in several studies |24ff27] . explaining recent nu- 
merical findings f26]. Our results show that the SL tran- 
sition at Af, for which we derived an explicit estimate 



(Eq. (Ill) is based upon two constituent parts. The first 
part is a ST crossover, which we studied both within 
and beyond the mean-field description. The second part 
is based on the probability, that the complex dynam- 
ics leads to a local energy above the PN energy barrier. 
Given the simplicity of the initial condition ^ used to 
probe the ST crossover, we expect that the crossover is 
not only experimentally readily accessible, but that its 
understanding could also be vital in generating long-time 
coherent states, without the need to fine-tune the initial 
state. 
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